On the faint end of the high redshift AGN luminosity function 
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ABSTRACT 

Using the results of recent optical surveys we conclude that the non-detection of quasars 
down to faint magnitudes implies a significant flattening of the high redshift (z ~ 6) optical 
active galactic nuclei (AGN) luminosity function for M1450 > —24.7. We find that all the data 
are consistent with a faint-end slope for the optical AGN luminosity function of (3 = —2.2 and 
j3 = —2.8, at the 90% and 99% confidence level respectively, flatter than the bright-end slope 
of f3' ~ —3.2. We also show that X-ray deep surveys have probed even fainter magnitudes than 
the optical ones yielding more significant constraints on the shallow faint-end slope of the optical 
luminosity function. The inclusion of Type II AGN candidates, detected in the Chandra deep 
fields, hints towards an higher normalization for the total AGN luminosity function, if these 
sources lie at 5 < z < 6.5. We then discuss simple theoretical models of AGN formation and 
evolution in the context of cold dark matter cosmology. The comparison with the total AGN 
luminosity function favors a redshift-dependent relation between black hole and dark matter halo 
masses of the type M, oc M^ alo , with 1.3 < a < 1.7, compatible with independent studies from 
statistical analysis and rotation curve measurements. Finally we compute the quasar contribution 
to reionization to be < 9% at z ~ 6, up to ~ 30% when integrated within 5.5 < z < 6.5, 
significantly smaller than that from galaxies. 



Subject headings: galaxies: active - quasars: general 

1. Introduction 

A fundamental tool to study the evolution of 
Active Galactic Nuclei (AGNs) is their luminosity 
function (LF) and its evolution with redshift. This 
enables us to probe the accretion history of Super- 
massive Black Holes (SMBHs) which, in turn, set 
constraints on galaxy evolution, as SMBHs are be- 
lieved to be ubiquitous and tightly linked with all 
stellar bulges and spheroids today (e.g. Ferrarese 
& Ford 2005 and references therein). Moreover 
constraining the AGN LF is essential to under- 
stand the AGN contribution to the total emissivity 
background in different bands, their environmen- 
tal impact, their clustering properties and the role 
AGNs had, if any, in the cosmological re-ionization 
of Hydrogen and Helium. 

X-ray observations carried out with Chandra , 
ASCA, HEAO-1 are statistically complete in X- 
ray to a level of ~ 90% up to z ~ 3 4 (e.g. La 



Franca et al. 2005). These observations have as- 
sessed that AGNs statistically follow a luminosity- 
dependent density evolution (e.g. Ueda et al. 
2003), in which lower luminous AGNs peak at 
z < 1, while higher luminous ones, usually op- 
tically bright, peak at higher redshifts, z ~ 2 (e.g. 
Osmer 2004). However how far in time can we 
track AGN evolution? The picture is still not yet 
complete at higher redshifts, where X-ray mea- 
surements of the AGN LF are lacking. 

On the optical side the SDSS team (Richard 
et al. 2006), and earlier works by Kennefick et 
al. (1995), COMBO-17 (Wolf et al. 2003) have 
probed the bright tail of the optical AGN LF up 
to z ~ 4.7 down to magnitudes of M1450 < —23.3 
showing that a simple power-law can easily fit the 
data. The GOODs team (Cristiani et al. 2004; 
Fontanot et al. 2006a) has been able to set con- 
straints on the shape of the optical LF down to 
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-^1450 ~ —21.5 and up to z ~ 5, showing evi- 
dence for a flattening at the faint-end. Fan et al. 
(2004; F04 from here on) have instead probed the 
very luminous tail of the z ~ 6 optical LF up to 
M 145 o < -26.5. 

Here we use existing observations of optically 
selected AGN and Chandra deep fields to deter- 
mine constraints on the redshift 6 shape of both 
the optical and X-ray LF. We will in particular 
show that X-ray observations allow us to probe 
the 2^6 LF to ~ 2 magnitudes fainter than the 
deepest optical observations. In section [2] we de- 
scribe our method, present the constraints on the 
faint end LF in section [3j their implications in sec- 
tions S] and [51 and conclude in section [5] We use 
concordance cosmological parameters £l m = 0.27, 
fl A = 0.73, Q h = 0.022/i" 2 , H = 71 km s" 1 
Mpc" 1 , a 8 = 0.84, and T = 0.2 (Spergel et al. 
2003). 

2. Method 

In the following we will always refer to the total 
population of active galactic nuclei as AGNs. We 
will call " QSOs" (or " quasars" ) the type I, optical 
AGNs, and Type II AGNs those which have a col- 
umn density above log-ZVy/cm -2 = 22. The latter 
kind of sources are often very faint or even missed 
in optical surveys, but still detected in X-rays. 

Several deep X-ray and optical surveys have 
now probed several portions of the AGN luminos- 
ity function at high redshifts, however often find- 
ing a scarce number of faint sources. In this paper 
we will show that the paucity of sources at these 
faint fluxes gives the most up-to-date stringent 
constraints on the knee and the faint-end slope 
of the QSO LF and give hints to the overall AGN 
distribution at z ~ 6. Similar studies, though with 
looser upper limits, have been carried out by Sharp 
et al. (2004), Willott et al. (2005), Wyithe (2004) 
and Dijkstra & Wyithe (2006). 

To convert X-ray to optical rest-frame luminosi- 
ties we will use the equation (see Richards et al. 
2005) 

logI/ 2 koV = logL2500 + Ckoxlog ( — ^ ) , (1) 

V ^2500 / 

between monochromatic luminosities and a con- 
stant photon index T — 1.8 to compute the broad- 
band X-ray luminosities. Many authors (e.g. Vig- 



nali et al. 2003) have found that the spectral index 
a ox depends significantly on luminosity. In par- 
ticular, the latest results by Steffen et al. (2006), 
who have collected a large number of optically and 
X-ray selected AGNs at several redshifts, claim a 
significant correlation of the spectral index with 
luminosity and a marginal one (l.lcr) with reshift, 
i.e. 

a ox = (-0.126 ± 0.013) logL 250 o - (0.01 ± (2) 
0.009)z+ (2.311 ± 0.372), 

where £2500 ' s ^ ne monochromatic luminosity at 
2500 A, in units of ergs -1 Hz -1 . 

In the following we will use results from sev- 
eral optical/NIR surveys carried out on different 
areas. To compare the luminosities probed by dif- 
ferent observations, we will first need to convert 
the broadband zab, where most of the observa- 
tions are carried out, to the apparent magnitude 
at wavelength 1450(1 + z)A, i.e. AB 1450 ( 1+Z ), and 
finally convert the latter to the rest-frame lumi- 
nosities Mi4 50 . 

We convert zab to AB 1450 ( 1+Z ) magnitude fol- 
lowing the spectrum described in Cristiani & Vio 
(1990) and updated in Cristiani et al. (2004). To 
convert then to absolute magnitudes M1450, we 
note that in the AB system ttiab = — 2.5 log /„ — 
48.6, where the monochromatic flux f v is in units 
of erg cm -2 s _1 Hz -1 , and the absolute magni- 
tude is Mab = — 2-51og_L„( 1+z ) + 51.6, where 
Ljy(i+2) is the luminosity at rest frame frequency 
v(\ + z), in units of ergs _1 Hz _1 . Using f v 
; 1 + z)L„( 1+z )/ '4:TtD^(z) leads to (see also Stern et 
al. 2000b) 

M 1450 = vLB 145 o(i+,) + 2 - 5 lo S(l + z) + 97.45 (3) 

-5 log D h (z) 

and (Fan et al. 2001) 

M B = M1450 + 2.5alog(4400/1450) + 0.12 (4) 

where Mb is the absolute B magnitude calibrated 
for a f v cx v a (where we set a = —0.5) spectrum 
and the factor 0.12 comes from the zero point dif- 
ference between the AB and Vega-based system 
(Fan et al. 2001). 

To constrain the shape of the z ~ 6 AGN LF, 
following Willott et al. (2005), we first define a 
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double power-law model LF of the type 

$(M U50 ,,) = ^%^ = p(z)x (5) 



dM 



$ 



lOOA((3+l)(M li50 -M* i5a ) + IQ0A(I3'+1)(M 1A50 -M^ 50 ) ' 

where p(z) = icr - 48 (^o), with z Q = 6.07, is 
the average redshift of the F04 LF. With such a 
parametrization we assume that the high-redshift 
AGN LF evolves with a Pure Density Evolution 
(PDE) beyond redshift of 5 at the same rate as the 
LF in Fan et al. (2001). Following the observa- 
tions by F04, $o nas been calibrated to match the 
cumulative number density of ~ 6.4 x 10~ 10 Mpc 3 
for sources with Mi4 50 < —26.7 at z — z (these 
are the limit and number density for Q m = 0.35, 
f^A = 0.65, Hq = 65 and have been converted 
to our cosmology). The bright-end slope and the 
break magnitude have been fixed to p — — 3 and 
M* 450 = —24.67 respectively, as this provides a 
very good match to the Cool et al. 2006 observa- 
tions (see below). The faint-end slope /3 has been 
left variable and constrained by matching the ob- 
served number density of sources in each survey. 
At z ~ 5 the LF in equation ([5]) is in agreement 
with the cumulative number density determined 
by Fan et al. (2001). 

Thus given a model LF we integrate it within 
the effective volume of the survey using the follow- 
ing equation (see e.g. Fan et al. 2001) 



7Jobs = Af2 



dz dV x 
dzdfl 



(6) 



/ $(L,z)dLp(L,z,SED) . 



and tune the faint-end slope in equation (JB|) to 
match the observed number of sources n Q b s in the 
survey. L m i n (z) refers to the redshift dependent 
absolute luminosity, corresponding to the mini- 
mum optical apparent magnitude (or the mini- 
mum flux) considered. We will in general assume 
a constant selection function p(L,z,SED) for each 
survey, as indicated by the authors or derived from 
their papers. In the following we will always use 
the form in equation ^ for <fr(L, z), unless in some 
specified cases, in which we will adopt the Fan et 
al. (2001) LF. The latter is a single power-law 
function valid in the range 4 < z < 5.5, in which 
we let the slope Pq to be variable, while we keep 



the redshift evolution and normalization at Mj" 450 
as in Fan et al. (2001; both M*^ 50 and normaliza- 
tion have been converted to our cosmology, which 
is required at each redshift bin during the evolu- 
tion). 

3. Results 

3.1. The Type I AGN Luminosity Func- 
tion 

A significant improvement with respect to pre- 
vious works has been to fix the knee of the LF 
in equation ^ to the value recently inferred by 
Cool et al. (2006). These authors discovered three 
objects down to zab=21.93, (i.e M1450 ~ —24.67 
at z= 5.53 using equation |4}) in the AGN and 
Galaxy Evolution Survey (AGES) spectroscopic 
observations of the NOAO Deep Wide-Field Sur- 
vey (NDWFS). They covered 7.7 deg 2 of the ND- 
WFS field with an average completeness level of 
54% down to Jab = 22 mag. These authors claim 
their observations agree with the F04 LF with a 
slope of p = —3.2 and a redshift evolution as in 
Fan et al. (2001) and in equation {6|). Using equa- 
tion ([7]), we find ~ 3 sources with zab ^ 21.93 
at z > 5, in very good agreement with observa- 
tions and with the authors. The integration was 
pursued up to z = 7, as beyond this limit the op- 
tical selection of quasars becomes impossible be- 
cause the Lya line redshifs out of the ZAB-band 
(see Cool et al. 2006). 

Following a similar procedure we find signif- 
icant evidence for a flattening of the QSO LF 
below A/1450 ~ -24.67. Recently, Mahabal et 
al. (2005, M05 hereafter) in the field around the 
z = 6.42 quasar SDSS Jl 148+5251 discovered a 
quasar at redshift z — 5.70 with luminosity of 
A/1450 ~ —23.6, i.e. zab = 23.0. This quasar 
is the faintest known optically selected quasar at 
z > 5.5, allowing determination of the faint end of 
the quasar LF to fainter luminosities. 

The authors claim the source, identified in a 
~ 2.5 deg 2 survey, has a surface density consis- 
tent with the Wolf et al. (2003) LF extrapo- 
lated through a PDE within their effective vol- 
ume, which follows the i?-band drop-out criterion 
of Stern et al. (2000), i.e. 4.3 < z < 5.8. We have 
checked that adopting the latter determination for 
the effective volume, we recover the results of the 
authors only if an average ~ 28% selection detec- 
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tion probability is used. However using the Wolf 
et al. (2003) LF, which is defined only for z < 4.7, 
needs an extrapolation to higher redshifts. More- 
over the functional form in equation ([B]) has been 
calibrated on the data by Cool et al. (2006) at 
z > 5. We therefore replace equation ([6]) with the 
Fan et al. (2001) QSO LF, more appropriate for 
the redshift bin probed by the i?-band drop-out 
technique. 

Keeping the same value for the selection func- 
tion we find that the detection of one object in 
the effective volume probed by M05, is consis- 
tent with a number density at M1450 ~ —23.6 of 

\og[dn/dM] 6.70 at the 99% confidence level 

(CL; with O = -2.8), down to \og[dn/dM] ~ 
-6.95 at the 90% CL (with O = -2.54; see Ta- 
ble [TJ we plot the M05 data in Figure [1] shifted by 
-0.15 magnitudes for avoiding overcrowding with 
Willott et al.'s data). Their data, shown as solid 
triangles in Figure ([1]) , are consistent with extrap- 
olating equation ([6]) to fainter optical magnitudes 

using = -2.3 and 1.7 at the 99% and 90% 

CL respectively (see Table [TJ. Note that these are 
conservative estimates. Extending the computa- 
tion of the LF down to the limiting magnitude of 

the M05 survey (z A b = 24.5, Afi 450 22) or 

using a higher completeness function would yield 
even lower limits on 0. 

Willott et al. (2005) undertook a multicolor 
quasar optical imaging with the Canada-France- 
Hawaii Telescope. They found a lack of quasars 
down to a limiting magnitude of zab =23.35 (i.e. 

M1450 23.63 at z = 6.1; see Figured]). The few 

detected quasar candidates had been subsequently 
rejected as being identified with low-mass stars in 
the near-infrared follow-up imaging. Using equa- 
tions ([6]) and ([7]), we find that the Willott et al. 
(2005) results are consistent with a faint-end slope 
of = -2.95 at 99% CL, down to = -2.27 at the 
90% CL. We use the selection function computed 
from simulated colors for high-redshift QSOs in 
the redshift range 5.5 < z < 6.7 (see Figure 2 in 
Willott et al. 2005). This is the most conservative 
selection function from Willott et al. (2005), de- 
termined for very extreme red color cuts. Using a 
less conservative estimate (solid line in their Fig- 
ure 2) our results drop to = -2.84 at 99% CL, 
down to = -2.2 at the 90% CL. 

The above findings are shown in Figure [TJ The 
gray area shows the data with its lcr-error bars by 



F04 for the z ~ 6.1 LF, while the solid dots show 
the extrapolation of the F04 LF to faint luminosi- 
ties assuming a single slope of = —3.2. The pres- 
ence of a significant downturn in the QSO LF for 
magnitudes below M1450 ~ —24.67 is evident. In 
particular the M05 observations are significantly 
stringent in constraining the faint-end slope. 

In Figure [TJ we also show the interesting result 
by Stern et al. (2000). These authors have discov- 
ered a quasar at z = 5.50 from deep, multicolor, 
ground-based observations covering 74 arcmin 2 
with zab = 23.4. To estimate their average se- 
lection function we redo the calculations in Stern 
et al. (2000) and integrate down to Mb — —22.5 
(in their cosmology) the Boyle et al. (1991) z = 2 
quasar LF, evolved as in Schmidt et al. (1995). 
The effective volume probed is again defined by 
the i?-band drop-out as in M05 (Stern et al. 2000). 
We find that an average completeness of p= 0.89 
recovers their quoted ~0.15 quasars in the field of 
view. As for M05, we then use the more up-to- 
date Fan et al. (2001) LF, more suitable for the 
redshift range probed by this survey, and equa- 
tion ([7]). Even in this case we let the slope to be 
variable, 0o = —2.92 in the specific case, while we 
keep the redshift evolution and normalization at 
Mf 450 as in Fan et al. (2001). 

We integrate down to ^4i?i45o(i+ z ) = 24.1 (see 
Table 1 in Stern et al. 2000), instead of using the 
measured zab = 23.4. This is because the au- 
thors find their source to have unusual spectral 
features for which they derive a special contin- 
uum model which is significantly different from 
the average spectral template used in this pa- 
per. The detection of one source in their field 
of view is consistent with a comoving density 

of \og[dn/dM (-22.5, 5.5)] 5.66 Mpc~ 3 mag" 1 , 

lower than the F04 extrapolation at the > 80% 
CL (see Figure [TJ . The number density we find 
still shows evidence for a turnover in the QSO LF 
expressed in equation ©, with = —2.8, though 
less significant than the above results, in agree- 
ment with the authors. 

We now turn to the analysis of X-ray data. 
Barger et al. (2003; see also Cowie et al. 2003) 
have used deep multicolor optical data to search 
for z > 5 AGNs in the 2 Ms X-ray Chandra 
Deep Field-North exposures. They found 500 
sources in their field of view out of which 423 have 
zab < 25.2. Color analysis was carried out on all 
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Table 1 

Faint-end slopes for the QSO LF at several confidence levels 



survey |, 


8\ 90% 


95% 


99% 


Area 


P 


Mahabal et al. 


1.7 


1.9 


2.3 


2.5 sq. deg. 


28 % 


Willott et al. 


2.27 


2.57 


2.95 


3.32 sq. deg. 


variable 


Barger et al. (single source) 


2.32 


2.38 


2.48 


160 arcmin 2 


50 % 



Note. — All the values in the first three columns give the faint-end slopes (3 
(equation @) needed to reproduce the observed number density in each survey; all 
confidence levels are computed through Poissonian statistics following Gehrels et al. 
(1986); the column before last indicates the effective area probed by each survey; 
the last column indicates the average selection function used for each survey. 



the sources, while 249 sources have spectroscopi- 
cally confirmed redshifts. Of these, only one ob- 
ject was found at 5 < z < 6 and none at z > 6. 
They found a faint AGN at z = 5.19 with zab = 
23.9, corresponding to an absolute magnitude of 
Mi45Q ~ —22.8, i.e. nearly a magnitude fainter 
than the M05 source, showing the power of X-ray 
observations. Barger et al. (2003) claim that none 
of the spectroscopically unidentified sources with 
zab < 25.2 is red enough to have z > 5. With the 
aid of the V/V max method the authors estimated 
the LF in the comoving volume 5 < z < 6.5 to be 
consistent with the extrapolation of the Fan et al. 
(2001) LF with a slope of (3 = -2.6 in the lumi- 
nosity bin 43 < logL 2 _iokcV < 44 (in ergs -1 ). 

Here we proceed with a similar calculation but 
using equation {2} and extrapolating the Fan et 
al. (2001) LF at redshifts beyond z ~ 5, in the 
same manner we computed the number density for 
Stern et al. (2000). We find that the number of 
observed sources down to zab = 25.2 integrated 
between 5 < z < 6.5, with j3 = —2.1 and a com- 
pleteness of p= 0.7 (Cowie et al. 2002) is ~ 1, in 
very good agreement with observations. We con- 
sidered an area of 160 arcmin 2 , corresponding to 
the solid angle where the hard X-ray flux limit has 
a signal-to-noise ratio of S/N=3 (see Figure 19 in 
Alexander et al. 2003). The X-ray flux consid- 
ered here is /o.5-8koV ~ 2.5 x 10~ 16 erg cm -2 s _1 
(Alexander et al. 2003; Barger et al. 2003). This 
flux limit is the minimum level of X-ray flux which 
corresponds to an optical magnitude brighter than 



zab = 26.9 (the optical detection limit in Barger 
et al. (2003)) within the whole effective volume 
considered. This is required as all the 500 sources 
in the CDFN have optical colors. 

We find a number density of log [dn/dM 
(-21.4,5.7)] ~ -5.70 (Mpc^mag- 1 ) correspond- 
ing to log L 2 -iokoV ~ 43.6 (ergs -1 ) (equation ©), 

a factor of 2 3 higher, but broadly consistent, 

with the results by Barger et al. (2003; plotted in 
Figure [T] as an open circle with its 99% CL error 
bars). 

Indeed deep X-ray observations effectively al- 
low us to probe even fainter optical luminosities. 
Mathur et al. (2002; Brandt et al. 2002) de- 
tected three z ~ 6 quasars discovered by Fan et 
al. (2001b) with Chandra with exposure times of 
less than 10 ksec. The objects are detected in 
X-rays up to 55 keV in rest frame, demonstrat- 
ing the incredible sensitivity of Chandra to de- 
tect faint sources. The faintest quasar in their 
sample, SDSS 1030+0524, at z = 6.28 was de- 
tected in about 8 ksec with an observed flux of 
/o.5-8kcV = 4 x 10 _15 ergcm _2 s _1 . The Chan- 
dra deep field observations have exposure time of 2 
Msec, so could detect down to 2000/8 = 250 times 
fainter sources with a flux limit of /o.5-8kcV ~ 
1.6 x 10 -17 erg cm -2 s _1 . This however would be 
the sensitivity limit at the aim point, close to the 
aim point sensitivity reported by Alexander et al. 
(2003). Therefore we expect the Chandra observa- 
tions to be sensitive in detecting sources down the 
minimum flux considered, i.e., we can extend the 
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computation using equations ([6]) and down to 
the limit of /o.5-8kcV ~ 2.5 x 10~ 16 erg cm -2 s _1 
(i.e. down to zab ~ 26 at the average redshift of 
z = 5. 7). In Figure Q] and Table [1] we report our 
findings. Using a completeness of about 50%, as 
computed by Cowie et al. (2002) for faint fluxes 
(notice that the use of an higher completeness level 
would strengthen our results), we find that the 
detection of a single source in the Chandra Deep 
Field in the range 5 < z < 6.5, is consistent with 
a faint slope of j3 = -2.48 at the 99% CL, down to 
/3 = —2.32 at the 90% CL (shown as solid squares). 

Barger et al. (2003) claim that none of the 
unidentified sources is red enough to be at z > 6. 
They however also point out that up to nine out 
of the optically faint spectroscopically unidentified 
sources with zab > 25.2, could still possibly lie in 
the range 5 < z < 6.5. Therefore assuming these 
sources to rely in the same luminosity bin as the 
single z = 5.19 source, they find an upper limit to 
their comoving LF about an order of magnitude 
higher than their previous estimate. Following the 
same approach undertaken for the case of a single 
source, we apply equations ^ and ([7]) down to 
jo.5-8kcV ~ 2.5 x 10~ 16 ergcm~ 2 s _1 and using a 
completeness of p = 0.5, proper for faint fluxes 
(Cowie et al. 2002). We consider all the nine 
sources not to be heavily obscured, i.e. they are 
Type I AGN. We find that the comoving number 
density of the LF to be log[dn/dA/(-20.6, 5.7)] ~ 
— 4.52 (Mpc~ 3 mag^ 1 ), using a slope (3 = —2.57 
(shown as an open square with its 99% CL in Fig- 
ured]). Our estimate is higher but still consistent 
with the one by Barger et al. (2003). To notice 
also that our estimate is also quite conservative, 
lower limits would be inferred extending the com- 
putation to lower magnitudes, close to the survey 
limit (z AB = 26.9). 

The 1-cr uncertainties in equation (TTJ) produce 
uncertainties of ~ 0.15 magn in the final compu- 
tation of the optical M1450 magnitude, as shown 
in Figure [TJ We also notice that our conclu- 
sions are quite insensitive to the actual fit used 
for the optical-to- X-ray flux ratios. For example, 
similar results are also obtained using the rela- 
tion between luminosity in the rest-frame R band, 
Lr, and the hard X-ray luminosity £2-10 keV cal- 
ibrated by La Franca et al. (2005; see their equa- 
tion (4)) for Type I AGNs. The Vignali et al. 
(2003) fit instead of equation (fl]) would yield even 



fainter optical magnitudes for the same X-ray flux 
and number density, thus strengthening our re- 
sults for a presence of a flattening at the faint-end 
of the QSO LF. 

We conclude this section noting that all the 
data are consistent with a break in the QSO LF 
at M!45o > —24.67. The faint-end slope of the op- 
tical Type I AGN LF, is consistent with being on 

average (3 2.2, at the 90% CL and (5 2.8, 

at the 99% CL (see Figure p. Our overall find- 
ings are consistent with the Fan et al. (2001) LF 
with (3 — —2.58 extrapolated to z = 6.1 (shown 
with "plus" signs in Figure [lj. This clearly shows 
how the high-redshift LF is more consistent with 
a PDE beyond z ~ 5. Richards et al. (2006) have 
recently found from a sample of SDSS QSOs, that 
the high-redshift QSO LF bright slope flattens at 
z > 2.5 from \/3'\ = 3.2 down to \/3'\ > 2.4 at 
z ~ 5. We suggest instead that a double power- 
law, steeper at bright magnitudes and shallower 
at fainter ones, is more indicative for the QSO LF 
at z ~ 6. Our results on the optical LF are in- 
stead in good agreement with those by Fontanot 
et al. (2006a), who compiled the QSO LF from a 
sample of GOODs and SDSS QSOs in the range 
3.5 < z < 5.2. They find significant evidence for 
a bright-end slope of \/3'\ = 3.3, a break in the LF 
at Mj" 450 ~ — 25 and a strong flattening below this 
limit. 

3.2. The Type I plus Type II AGN Lumi- 
nosity Function 

In the previous section we have considered the 
detected sources in the CDFN as Type I AGNs. 
The single high-redshift source detected with spec- 
troscopic redshift of z = 5.19 can in fact be 
safely considered as a Type I AGN. Converting 
its measured hard X-ray flux /2-8kcV ~ 7 x 
10~ 16 erg cm -2 s _1 to optical magnitudes (using 
equation we get zab ~ 23.6, very close to its 
actual magnitude of zab ~ 23.9. However some of 
the nine optically faint spectroscopically uniden- 
tified objects found by Barger et al. (2003) could 
be considered moderately obscured sources, with 
very faint optical fluxes (though still detectable). 
In this section we give constraints on the faint-end 
shape for the AGN LF, including the CDFN spec- 
troscopically unidentified sources assuming that a 
certain fraction of them are Type II AGNs lying 
at z ~ 6. 
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Indeed the 2 — 8keV band corresponds to the 
14 — 56 keV intrinsic band at z = 6, which is 
not much affected by absorption for column densi- 
ties log Nh < 24 (cm -2 ). Here we notice however 
that we do not expect the detection of any very 
obscured, Compton-thick source in the CDFN, 
with logL 2 -iokcV ^5 44 (ergs -1 ). From Wilman 
& Fabian (1999) we have taken the spectrum of 
two typical AGNs with column densities log Njj = 
24.25 (cm -2 ) and \ogN H = 24.75 (cm -2 ), derived 
from a Monte-Carlo realization which includes 
Compton down-scattering which reduces the in- 
tensity even in the highest bins of energy. We 
then computed the ratio between absorbed and 
unabsorbed luminosities (the former is the actual 
observed flux) for the various spectra in the z = 6 
redshifted, restframe energy hard band 14—70 keV 
(the least absorbed band of the CDFN). 

We find that the ratio L^fL 70kcV /L^ b 7 s 0kcV is 
- 71% for \ogN H = 24.25 (cm -2 ) and ~ 30% for 
log Njj — 24.75 (cm -2 ) (to notice for comparison 
that an AGN with logN H = 23.5 (cm -2 ) has a 
ratio of ~ 96%). A typical source with intrinsic 
luminosity L2-iokcV = 10 44 (ergs -1 ) (fainter lu- 
minosities would be below the X-ray flux limit) 
in the observed frame will in turn produce ob- 
served fluxes /2-8koV ~ 2 x 10 -16 erg cm -2 s -1 
and /2-8kcV ~ 1 X 10 -16 erg cm -2 s -1 for respec- 
tively the log N H = 24.25 (cm -2 ) and logN H = 
24.75 (cm -2 ) case. Therefore sources with log Nh 
> 24 (cm -2 ) and logL 2 _iokcV < 44 (ergs -1 ), 
would be hardly visible in the CDFN field of 
view (we remind that the flux limit is /2-8kcV ~ 
2 x 10 -16 ergcm -2 s -1 ). 

On the other hand moderately obscured sources 
(with log Njj ;$ 24 (cm -2 )) would definitely be 
visible at these high redshifts, and some of them 
might have been already discovered in the deep 
Chandra fields, among the optically faint Barger 
et al. (2003) sample. Other possible indications of 
the presence of Type II sources in the high redshift 
universe also come from Koekemoer et al. (2004). 
From the 2 msec Hubble Deep Field-North and 
the 1 msec CDFS, these authors have robustly 
detected a sample of seven sources with extreme 
optical/X-ray flux ratios >10 at a flux limit above 
/o.5-8kcV ~ 3 x 10 -16 ergcm -2 s -1 . Although 
these kind of sources are generally rare, if they 
are AGNs lying at z > 6 such that their Lya 
emission has been redshifted out of the zab band- 



pass, then they would definitely be high-redshift 
obscured sources. We cannot directly make a pre- 
diction with their data as we are converting to 
optical luminosities and these are special sources, 
strong outliers with respect to equation §5§ (also a 
conversion to bolometric luminosity would be dif- 
ficult in this case). 

In this section we still adopt the same X-ray K- 
correction as for unabsorbed sources with T = 1.8 
to convert flux to luminosities. We still use equa- 
tion ([3]), which holds for Type I AGNs, for con- 
verting X-ray fluxes to intrinsic optical luminosi- 
ties for Type II AGNs, assuming that unifica- 
tion models are on average valid. First we con- 
sider all the sources detected by Barger et al. 
(2003) as optically faint Type II AGNs lying at 
z > 5. We then integrate equation ([6]) in the 
redshift range 5 < z < 6.5 above a flux limit 
of /2-8kcV — 2 x 10 -16 erg cm -2 s -1 correspond- 
ing to an effective area of 40 arcmin 2 (see Figure 
19, bottom panel in Alexander et al. 2003). The 
adopted flux limit corresponds to apparent mag- 
nitudes zab > 25.5 at z ~ 5.7 (equation Q). Us- 
ing f3 = 3.45 we obtain ~ 9 sources in the field 
of view, corresponding to a comoving density of 
log[dn/dM(-21.2,5.7)] ~ -3.85 Mpc -3 mag" 1 . 

The results are shown with an open square with 
its 99% confidence level in Figure [2] with its error 
bar on the optical luminosity when using equa- 
tion ([3]). If the nine sources are confirmed to be 
Type II AGNs at z ~ 6, this would show evidence 
for an extraordinary contribution of Type II AGNs 
at the faint end of the AGN LF. Assuming all these 
faint sources to lie at z ~ 6 would imply a fraction 
obscured-to-unobscured AGNs of order 9. X-ray 
surveys at lower redshifts found much more mod- 
erate numbers for this ratio (e.g. Ueda et al. 2003, 
La Franca et al. 2005, Gilli et al. 2006). 

In fact the above scenario is an extreme one, 
as most probably many of these optically faint 
sources probably belong to lower redshifts (see dis- 
cussions in Barger et al. 2003 and Cowie et al. 
2003). Therefore we more safely assume that up 
to 3 sources, out of the 9 detected, i.e. three times 
the surface density of the Type I LF (made up 
of only the single z = 5.19 source), are actual 
contributors to z ~ 6 AGN LF with logN H < 
24cm -2 , following the results by La Franca et 
al. (2005), who find an almost flat distribution 
of the fraction of sources per logarithm bin of Njj . 
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In this second case we get a comoving density 

of log[efn/dM(-21.2,5.7)] 4.3Mpc -3 mag- 1 , 

shown as a solid square in Figure [2] with its 99% 
CL error bars. 

We then compare our AGN LF estimate with 
what would be obtained by multiplying the av- 
erage QSO LF (with a faint slope of /? = —2.5; 
section l3~Tj) by a factor 1 + R, where R is the ratio 
between Type II and Type I AGNs, as calibrated 
by several groups for AGNs detected at lower red- 
shifts. For computing R we also include sources 
up to log Nh < 26 cm -2 . We use the results by 
Ueda et al. (2003), La Franca et al. (2005) and 
Gilli et al. (2006). As specified above, we consider 
Type II AGNs all those sources with log Nh > 
22 cm -2 . Therefore by definition the number ra- 
tio Type II-to-Type I will be given by R = n[22 < 
log(7V ff /cm -2 ) < 26]/n[log(^/cm- 2 ) < 22]. 
For Ueda et al. (2003) and La Franca et al. (2005) 
we use their quoted distributions for the fraction 
of objects up to log Nh — 26 (cm -2 ). 

Gilli et al. (2006) developed an updated model 
for the synthesis of the X-ray background, for 
which they adopted a different definition for the 
number ratio between obscured and unobscured 
AGNs, i.e. R = »[21 < log(^/cm -2 ) < 
24]/n[log(A r ff/cm -2 ) < 21]. According to their 
model this ratio is redshift-independent, luminosity- 
dependent and, translated to optical magnitudes, 
it reads as 

R(M 1450 ) = R s x exp (_io -0 - 4 (^«o - ^)) + (7) 

R q x (l-cxp(-10 -0 - 4 ( M " 50-M ^)) , 

with R s = 3.6, R q = 1 and M c = -20.7, 
i.e. log£ .5-2kcv ~ 43.5 (ergs -1 ), where R s 
and R q are the n[21 < log{N H /cm -2 ) < 24] 
/n[log(iV^f/cm -2 ) < 21] ratio at the low- and 
high-luminosity regimes respectively, and M c = 
—20.7 is the characteristic luminosity dividing the 
two regimes (Gilli et al. 2006). To include also 
Compton-thick sources with column densities up 
to logNn — 26 (cm -2 ), we consider double this 
ratio (see Figure 7 in Gilli et al. 2006). 

The resulting AGN LF from the Gilli ct al. 
(2006) model is shown with a dashed line in Fig- 
ure [2] In the same figure we also plot the results 
of using R as calibrated from Ueda et al. (2003; 
dot-dashed line) and La Franca et al. (2005; long- 
dashed line). All the models for the distribution 



of sources as a function of Nh considered provide 
similar estimates for the total AGN LF, and are 
consistent with the extrapolated data at the faint 
end, when using 3 sources as proxy for the z ~ 6 
AGN population within the effective area. We re- 
mind here that our estimate for the faint-end AGN 
LF could easily be underestimated, as the CDFN 
cannot detect Compton-thick sources (those with 
log N H > 24cm -2 ). 

To notice that all the models also provide a sig- 
nificant number of sources at the bright end of the 
AGN LF in excess of the QSO LF. Due to the small 
volume probed by the CDFN we however cannot 
provide any kind of constraint on the actual con- 
tribution of intrinsically very luminous and more 
rare obscured sources. Anyway in the following 
we will consider as our best estimate of the to- 
tal AGN LF at z ~ 6, the QSO LF multiplied by 
(l+i?Uoda), where i?uoda is the Type II-to-Type 
I ratio as extracted from Ueda et al. (2003). In 
the following section we will compare theoretical 
models with such a determination of the AGN LF. 

4. Implications 

The high redshift AGN LF is a very sensitive 
tool to probe the connections between supermas- 
sive black holes (SMBHs) and their host dark mat- 
ter (DM) halos. In this section we discuss possi- 
ble physical implications for SMBH evolution from 
our estimated AGN LF. We compute the total 
AGN LF from the formation rates of DM halos 
simply as (e.g., Mahmood et al. 2005) 

x,it \ 4. r> r*r \ dM halo dM. 

<S>{L B , Z) ~ t A GN X i?MM(Mialo, Z) X — — — . 

dM, aL b 

(8) 

Here i?Mivi(-Mhaio, z) is the formation rate of DM 
halos with mass A/haio at redshift z, massive and 
virialized enough to host a galaxy /AGN; £agn 
in this model is the average timescale during 
which an AGN is active (but not necessarily visi- 
ble). We neglect any dependence on mass for this 
timescale and just set it equal to £agn = 10 7 yr. 
M, is the SMBH mass producing the luminosity 
Lb- In the following for simplicity we will use 
extended Press-Schechter (EPS) theory to com- 
pute i?MM(-Mhaio, z) and we will consider only ha- 
los logMhaio/Af© < 13.2. The limiting mass of 
Mhaio ~ 1O 13 M0 as host for the most massive 
galaxies and most luminous quasars has been con- 
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strained through clustering (e.g. Croom et al. 
2005), statistical (e.g. Shankar et al. 2006), and 
semi- analytical (e.g. Granato et al. 2004) tech- 
niques. 

Most probably the EPS theory is a rather poor 
approximation for describing the statistics of high 
redshifts AGN host halos, moreover a more physi- 
cal treatment of the baryonic processes is required 
when comparing with observations (e.g. Lapi et 
al. 2006, Fontanot et al. 2006b). However here 
we neglect such issues and focus on the most suit- 
able relation M, — Mhaio, among those suggested 
in the literature, which better compares with ob- 
servations within the very basic model expressed 
in equation ([5]). We therefore consider different 
recipes for the jacobian relation in equation ([8]) 
relating DM and SMBH mass. 

We first discuss the Wyithe & Loeb (2003; 
see also Mahmood et al. 2005) model which 
connects the SMBH and DM halo masses as 
M. cx M^f (l + zf/ 2 , derived from a basic AGN 
feedback-constrained model. We shall refer to this 
model as model A in the rest of the paper. Using 
Model A, we calculate the z = 6.1 LF, shown as 
a solid line in Figure [3] This model is consistent 
with the AGN LF estimated in section 13.21 

Several other groups have also formulated the- 
oretical models for describing SMBH evolution in 
the context of cosmological theories (e.g. Haenelt 
& Rees 1993, Cavaliere & Vittorini 2000, Vittorini 
et al. 2005, Hopkins et al. 2005). Here we adopt 
a simple version of one of these models, viz. the 
model by Kauffmann & Haenelt (2000; see also 
Croton 2005) who parameterize the link between 
SMBH mass and halo mass as: 



M, = 0.03m R 



1 + 



/28()kit)s 



n^cold • 



(9) 

where wir represents the merger mass ratio be- 
tween the central and satellite galaxies, i.e. m R = 
TOsat/wcentrai, Kdr is the virial velocity of the AGN 
host halo and m co id is the amount of cold baryonic 
matter (which will eventually form stars). The 
equation ([9|) prescription with rriR = 1, consisting 
of only major mergers, is our model B. 

We then consider three cases in which the cold 
fraction is m co id = 0.1 x Afhaio, cricoid = 0.016 x 
Afhaio, and m co id = 0.005 x Afi ia i . The result- 
ing LFs are plotted in Figure [4] The first case, 



in which the amount of baryons per halo is close 
to the cosmological value, lies well above the ob- 
served LF. The second case fits the bright-end LF, 
assuming a baryonic fraction of a few percent, con- 
sistent with results by Shankar et al. (2006). It is 
interesting to note that a strong decrease of bary- 
onic mass fraction with decreasing halo mass is 
needed to reproduce the AGN LF, down to only 
0.5% at the very faint end. 

We now test a model based on Ferrarese et al. 
(2002; see also Shankar et al. 2006). We call this 
model C, and is the most empirical way to track 
the Af. ~ Afhaio relation. We relate the SMBH 
mass to bulge velocity dispersion using the Af. — a 
relation. The a is then connected to the circular 
velocity of the host galaxy by using the empirical 
a — Vc relation (from Baes et al. 2003). The cir- 
cular velocity is then scaled to the virial velocity 
assuming a Navarro, Frenk & White (1996; NFW) 
profile and a concentration parameter c{M\ ia \ 0l z) 
from Bullock et al. (2001). Finally the virial ve- 
locity is converted to halo mass using the virial 
theorem, GAfh a io/-Rvir = K 2 ir( z vir) (we use the re- 
lations in Loeb & Peebles 2003 given for each z V i r ), 
which is computed at the moment of the virializa- 
tion of the host DM halo (i.e. we set z v ir = 6.1). 
We also add a scatter of ~ 0.35 in the resulting 
Af. — Mhaio relation, as empirically suggested by 
Ferrarese (2002) and theoretically supported by 
Lapi et al. (2006). Our result reads as 



Af. = 5.8 x 10 7 



Afi 



halo 



10 12 Mi 



1.39 



1 



7 



Af 



© ) 



(10) 

which is close to the one calibrated locally through 
statistical arguments by Shankar et al. (2006) and 
to the one obtained by Granato et al. (2004) and 
Lapi et al. (2006) in semi-analytical models. 

One uncertainty in the above computations is 
the radius at which one measures the circular ve- 
locity. Ferrarese (2002) uses R ~ R s (also used 
in equation [TUJ) , where R s is the scale radius of 
the NFW DM profile. In Figure[5]we compare the 
predictions of model C with the calibrated high- 
z AGN LF. The model matches the data when 
R ~ (0.5 l)R s and it has the advantage of be- 
ing empirically-based, the only assumption here is 
that the M. — a and a — V c relations are almost 
constant in time. 

Our findings therefore differ from those in 
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Willott et al. (2005b), who attempted a simi- 
lar calculation and found that the number of DM 
halos at z ~ 6 falls short by several orders of 
magnitude in accounting for the number of high- 
redshift AGNs. Their conclusions were based on 
the use of the virial theorem calibrated at v ; r = 0, 
which inevitably leads to significantly underesti- 
mating the BH mass (and therefore AGN luminos- 
ity) for a given halo. Our findings are consistent 

with a ~ (0.5 0.6)V vil (M halo , z). The latter 

relation in turn has been used by Cirasuolo et al. 
(2005) to show that an impressive fit to the dis- 
persion velocity function (Sheth et et al. 2003; 
Shankar et al. 2004) is recovered when integrat- 
ing the DM formation rates Rmm (-^haloi ^vir) and 
setting a ~ 0.57 x K,ir(^haio, Aar) (see also Loeb 
& Peebles 2003 and Shankar et al. 2006). 

The simple approach undertaken in this section 
shows that a, redshift-dependent, super- linear re- 
lation of the type M. oc M^ alo , with 1.3 < a < 1.7, 
as the ones proposed in model A and C, is best 
suited to represent the AGN LF determined in 
section 13.21 Alternative scenarios would require 
a very strong variation of baryonic content and/or 
AGN activity timescale with halo mass. 

5. Contribution to Reionization 

In this section we compute the possible con- 
tribution of AGNs and stars to the reionization 
of the universe. In the following we will consider 
only Type I AGNs as contributors of the reioniza- 
tion along our line of sight. The transition from a 
neutral to a fully ionized Inter Galactic Medium 
(IGM) is usually statistically described by a differ- 
ential equation for the time evolution of the vol- 
ume filling factor of the medium Q(z). The latter 
quantity quantifies the level of the IGM porosity 
created by the HII, Hell and Helll ionization re- 
gions around radiative sources such as QSOs and 
galaxies. 

When the clumping factor of the medium is 
high, C » 1, which is always a good approx- 
imation at high redshifts (the simulations by 
Gnedin & Ostriker 1997 give C ~ 30), follow- 
ing Madau et al. (1999) one can approximate 
Q{t) ~ hi on /nHtrec, where t rcc is the volume- 
averaged mean recombination rate, rii on is the 
total rate of ionizing photons and fiu is the mean 
hydrogen density of the expanding IGM. 



The condition for reionization will then read 
Q w 1, which translates into a total required ion- 
izing photon rate (Madau et al. 1999; Fan et al. 
2001): 

(ii) 

We compare such a rate with those from QSOs 
and galaxies at similar redshifts. We start comput- 
ing the photoionization rate from QSOs, which is 
given by (e.g. Madau et al. 1999, Stiavelli et al. 
2004) 

n-AGN = / t dv , (12) 

Juh h P V 

where hp is the Planck's constant, vg — 3.2 x 
10 15 Hz is the frequency at the Lyman limit (912 
A, i.e. one Rydberg). The volume emissivity of 
the total population of sources is given as 

/oo 
<P{L,z)L u {L)dL, (13) 

where L v is in ergs Hz . The integration in 
equation (fT2|) is usually performed up to 4 Ry- 
dbergs, as photons more energetic than that are 
preferentially absorbed by He atoms. Anyway we 
have checked that extending the computation to 
~ 40 Rydbergs for example, increases our final 
results by just ~ 12%. L m - ln in equation (TT3]) is in- 
stead the minimum luminosity probed by the sur- 
vey, which we fix to M 1450 = —20.6, following the 
results in section |3~T1 

We thus find that the AGN contribution to 
reionization at redshift z ~ 6 integrating down 
to M 1450 = —20.6, varies from 2% to ~ 9%, us- 
ing a faint-end slope f3 of -2.2, -2.8 respectively. 
These results are rather independent of the de- 
tails for the AGN spectrum. For the computa- 
tions above we have used the spectrum described 
in Schirber & Bullock (2003), however using, for 
example, the three-power law spectrum adopted 
by Madau et al. (1999), reduces the AGN contri- 
bution by < 13%. 

We follow Meiksin (2005) for computing the 
photoionization rate from galaxies 

, dN phot , . 

™Galx — /csc-TTT Pstar ■ (,-1-4) 

"-'"star 

Here dN v ^ ot /dM star w 10 61 phot s _1 is the produc- 
tion rate of ionizing photons per solar mass of stars 
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formed, p sta r is the Star Formation Rate (SFR) in 
units of M Q yr _1 Mpc -3 , and / csc is the average 
fraction of ionizing photons which can freely es- 
cape the star-forming regions into the IGM. 

The photons per solar mass has been taken from 
Smith et al. (2002; see also Meiksin 2005) for a 
Salpeter Initial Mass Function (IMF) and a metal- 
licity of 20% the solar. The SFR was instead re- 
cently derived by Richard et al. (2006) to be ~ 
0.027M Q yr^ 1 Mpc~ 3 from a deep survey of lens- 
ing clusters aimed at constraining the abundance 

of star-forming galaxies at redshift z ~ 6 10. 

We then take results from the recent paper by 
Shapley et al. (2006), who have constrained the 
escape fraction in the ionizing Lyman continuum 
from a sample of 14 z ~ 3 star-forming galaxies, 

finding an average / csc ~ (11 18)%, and we 

use the average value of / csc ~ 15%. In the fol- 
lowing we neglect the diffuse emission from the 
IGM itself which however could add an additional 
photoionization rate (up to 50%, e.g. Meiksin & 
Madau 1993) over the direct one from the sources. 
We find that galaxies can contribute about 50% at 
z = 6.1 to reionization. 

We now take a step forward and integrate the 
contribution to reionization from both type of 
sources within the range of the Chandra deep 
fields, i.e. 5.5 < z < 6.5, and compare the in- 
ferred number of ionizing photons with the num- 
ber of available neutral hydrogen atoms at z ~ 6, 
n H = 6.2 x 10 66 (see Stiavelli et al. 2004). We 
find that the contribution from AGNs to grow up 

to (7 30)% while the one from galaxies being 

about 70%, keeping the SFR constant in this red- 
shift bin, as suggested by observations. The PDE 
evolution of the AGN LF for z < 6 makes the 
contribution from AGNs grow more than the one 
from galaxies. 

The photons produced by galaxies reach an av- 
erage contribution to reionization of > 70%, im- 
plying a minor contribution from QSOs, a conclu- 
sion in very good agreement with previous calcu- 
lations (e.g. Meiksin 2005, Fan et al. 2001, Barger 
et al. 2003). However we also notice that Bouwens 
et al. (2005) from the deep NICMOS fields in the 
Hubble Ultra Deep Fields (HUDF) have detected 
5 objects in 5.7 sq. deg., inferring a SFR about a 

factor of 2 3 below the one found by Richard 

et al. (2006). This would drop the galaxy contri- 
bution down to (25 40)%. Similar points were 



also recently addressed by Mannucci et al. (2006) 
who claim that the high-redshift SFR could fall 
short of completely reionizing the universe. They 
however also point out that the Bouwens et al. 
(2005) estimate should be further increased by a 
factor of 2 for dust extinction and/or accounting 
for luminosity evolution in the estimate. Using 
therefore the Bouwens et al. (2005) result in equa- 
tion (fT4|) would still leave room for a significant 
contribution from QSOs. Other ways to fill the 
possible " missing" photons from galaxies could be 
to alter significantly the photon luminosity from 
massive stars, invoking extra IGM diffuse light, 
or including other photoionizing sources (such as 
Population III stars and/or miniquasars powered 
by intermediate-mass BHs). However the contri- 
bution from the latter sources must be < 30% at 
the most. 

6. Conclusions 

In a recent paper, Cool et al. . (2006) found 
three QSOs at high redshifts (z > 5.4) down to 
magnitudes of M1450 ~ —24.7. The detection of 
these quasars is roughly consistent with the ex- 
trapolation of the bright steep slope of the quasar 
LF by F04 with (J = —3. At lower magnitudes not 
many QSOs are observed, as recently concluded by 
the optical/NIR surveys conducted by M05, Stern 
et al. (2005) and Willott et al. (2005). In this pa- 
per we show that such now-detection of QSOs im- 
plies a significant downturn in the optical QSO LF 
at lower magnitudes with respect to those probed 
by Cool et al. (2006). Moreover we show that 
the Chandra deep field observations probe about 

1 2 magnitudes deeper than the optical ones 

(down to M1450 = —20.6). The inferred faint end 
co-moving number density is lower by at least an 
order of magnitude with respect to the extrapola- 
tion of the F04 LF, at a high significance level. 

If a fraction of the nine Chandra deep field 
sources are Type II AGNs detected at X-ray fluxes 
/2-8koV ^ 2 x 10 -16 ergcm~ 2 s _1 and with very 
faint optical magnitudes (zab > 25.2), this would 
place a constraint on the shape of the Type I plus 
Type II AGNs at faint magnitudes, being consis- 
tent with a slope close to the extrapolation of the 
F04 LF, i.e. \0\ > 3. This estimate of the AGN LF 
is however representative of only the sources with 
\ogN H ~ 24 (cm" 2 ), as more obscured sources 
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would be below the CDFN flux limit. X-ray sur- 
veys at lower redshifts (e.g. La Franca et al. 2005) 
have shown that an almost flat distribution in the 
number of sources is observed up to column den- 
sities of log Nh ~ 24 (cm -2 ). We can therefore 
safely assume that only 3, out of the nine optically 
faint sources to be AGNs at z ~ 6. This yields a 
total AGN LF consistent with updated synthesis 
models for the XRBG. 

We then compare this result with predictions 
from three simple physically and empirically mo- 
tivated basic models for AGN evolution in DM ha- 
los. We generally find that a redshift-dependcnt 
relation of the kind M. - M£ alo , with 1.3 < a < 
1.7, produces models more consistent with the 
AGN LF. Such relations were also found from 
independent statistical studies by Shankar et al. 
(2006) and rotation curve extrapolations in spi- 
rals and ellipticals (Ferrarese 2002). Alternative 
models, such as varying the amount of cold gas 
fraction or AGN active timescale in smaller ha- 
los, would require extreme fine tuning and a large 
variation in the parameters with dark matter halo 
mass. Here we also notice that if we reduce the 
AGN duty-cycle just to the visible optical phase, 
i.e. substituting £agn with £qso i n equation ((5]), 
the match between models and the Type I AGN 
LF derived in section 13. 1\ would yield similar re- 
sults on the type of M, — Mhaio relation. 

Finally we have shown that AGNs can account 
for at most 30% of the total reionization of the 
Universe at redshifts around z ~ 6. However we 
also notice that a strong caveat on the actual role 
of galaxies in reionizing the universe, is set by the 
normalization of the cosmological SFR at high red- 
shifts, which is still uncertain by factors of about 
2. 
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Fig. 1. — Type I AGN Luminosity Function at z ~ 6. The squares represent the estimated AGN density at the optical 
magnitude of M1450 ~ —20.6 needed for a single source to be detected above the broad band X-ray flux limit of /0.5— 8kcV ~ 
2.5 X 10 — 16 erg cm -2 s — 1 in 160 arcmin 2 . From top to bottom the two solid squares are the 99% and 90% confidence level 
estimates, while the open square has been estimated taking into account all the nine spectroscopically unidentified sources in 
the deep Chandra fields with its 99% error bars; the horizontal error bar reflects the uncertainty in the computation of the 
optical magnitude at fixed X-ray flux when the l-cr uncertainty in equation (|3j are taken into account. The open circle is the 
AGN luminosity function estimated from Barger et al. (2003) with its 99% CL. The solid pentagon is the AGN density from 
Stern et al. (2000) with the 99% confidence level error bars; the solid diamonds are the AGN densities from Willott et al. 
(2005) again shown with 99%, 95% and 90% confidence level (from top to bottom); the solid triangles are the AGN densities 
computed from Mahabal et al. (2005), with 99%, 95% and 90% confidence level (from top to bottom); the open diamond is the 
AGN density computed from Cool et al. (2006) with its 99% confidence level error bars; the crosses are the Fan et al. (2001) 
luminosity function extrapolated to z = 6.1 assuming PDE from 2 = 5. The values reported on the upper-left of the figure are 
the absolute values for the faint slopes of the model AGN luminosity function in equation J6j; the shaded area is the observed 
luminosity function by Fan et al. (2004) with its 1 — a uncertainty while the thick-dotted line is its extrapolation to fainter 
magnitudes. 
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Fig. 2. — Type I plus Type II AGN luminosity function. The open square is the AGN luminosity function considering all the 
nine sources in the Chandra deep fields as Type II AGNs; the solid square is instead the AGN luminosity function computed 
when 3 sources are taken from the Bargeret al. (2003) optically faint sample to be z ~ 6 AGN candidates. The triple dot-dashed 
line is the Type I AGN luminosity function with the average slope of (3 = —2.5 and the striped area is the range of 90—99% CL 
for the Type I AGN luminosity function; the dashed line is the Type I AGN luminosity function multiplied by 1 + 2R where R 
is the Type II/Type I ratio as given in Gilli et al. (2006; equation J5)), while in the dot-dashed and long-dashed lines R has 
been extracted from Ueda et al. (2003) and La Franca et al. (2005) respectively. All the other symbols and lines are as in 
Figure [JJ 
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Fig. 3. — Estimated AGN luminosity function at z ~ 6 obtained by considering the average QSO luminosity function derived 
in section f3 , 1 1 ( with its uncertainty region shown as a striped region) and multiplying it by (1+-Rucda)> to account for all sources 
up to log Njj = 26 (cm -2 ). The solid line is the result of Model A (see text); all sources are assumed to radiate at the Eddington 
limit. 
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Fig. 4. — The dot-dashed, dashed and soiiei lines are the results of Model B (see text) for different fraction of cold gas with 
respect to the dark matter host mass as labeled; all sources are assumed to radiate at the Eddington limit. All the other curves 
are as in Figure [3] 
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Fig. 5. — The dot-dashed, dashed and solid lines are the results of Model C (see text) for different R/Rs ratios as labeled, 
assuming that all sources are radiating at the Eddington limit. All the other curves are as in Figure [3] 
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